in 
in 



Symmetric Triple Well with N on- Equivalent 
Vacua: Instantonic Approach 

Q : H. A. Alhendi^ and E. I. Lashin^'^'^ 

Csi ; ^ The Abdus Salam ICTP, P.O. Box 586, 34100 Trieste, Italy 

Q_i' ^ Department of Physics and Astronomy, College of Science, 

.^ ■ King Saud University, Riyadh, Saudi Arabia 

^ Department of Physics, Faculty of Science, 
Ain Shams University, Cairo, Egypt 

> ■ February 1, 2008 

O' 

O 

(N 

\o _ 

c7^ . We show that for the triple well potential with non-equivalent vacua, 

^T ' instantons generate for the low lying energy states a singlet and a doublet 

^^ \ of states rather than a triplet of equal energy spacing. Our energy split- 

(-H ■ ting formulae are also confirmed numerically. This splitting property is 

O ' due to the presence of non-equivalent vacua. A comment on its generality 

to multi-well is presented. 
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Abstract 



Instantons are non-trivial classical solutions of Euclidean field equations for 
^^ , which the action is finite Pj. Their importance, besides being topological con- 

$H ' figurations, conies from their finite contributions to Feynmann path integral. In 

quantum mechanics instantons correspond to non-trivial finite classical solutions 
of classical equations of motion with inverted potential |2] . 

The uses of instanton calculations have proven to be useful in analyzing non- 
perturbative aspects of quantum mechanical systems with degenerate vacua, this 
is because instanton solutions contribute to the quantum tunnelling phenom- 
ena, which can be calculated with the aid of the dilute gas approximation. A 
celebrated example is the splitting of energy level in the symmetric double well 
case 12 ■ 

To our knowledge only a few recent papers attempted to apply the instanton 
method to the triple well potential with non-equivalent vacua |3101|S]. This 
problem is rather involved compared to the symmetric double well case. It 
has been suspected that in the presence of non-equivalent vacua the dilute gas 
approximation may break down |S]. Moreover it has been claimed |31 0] that 
the average of the harmonic frequencies over the non-equivalent vacua of the 
potential serves as the central position for the equidistant nearly degenerate 



first three levels. Unfortunately this claim does not lead to the correct spectra 
and is in contradiction with the numerical calculations. Thus understanding 
the structure of the energy levels in a triple well becomes important. These 
issues encourage us to look at the problem more carefully and to our knowledge 
a proper treatment has not been previously carried out. 

In this work, we consider the triple well potential of the form, which admits 
inversion symmetry (x — > —x), 
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The potential ,V{x), shown in Fig. Q] consists of central, left and right well 
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Figure 1: Triple well: V(x) ■■ 
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separated by barriers. In the vicinity of each well (near the minima at x = 0, ±1) 
the potential can be approximated as a harmonic oscillator potential. The 
frequency corresponding to the central well is wq = ^i while for the left and 
right well is wi = 2 a; as can be easily verified by expanding the potential 
around each minima. These different frequencies show the different curvatures 
at each minima. 

Computation of transition probability amplitudes between different minima 
allows to extract the low lying energy eigenvalues. Here the instanton contri- 
butions to these transitions are calculated by making a proper treatment for 
the non-equivalent vacua. As will be shown this proper treatment has a drastic 
effect on the energy spectra in contrast with what has been claimed. 

For the potential in Eq. (Q we have four kinds of instanton solutions con- 
necting neighboring minima and have the forms: 
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where /(/) indicates instanton(anti-instanton) respectively. 



The transition from the minimum at a; = to that at a; = 1 during the 
EucHdean time interval T can be written as: 

(l|e-^^|0) = {l\E,) {E,\0) e-^°^ + (l|i?2) (i?2|0) e'^-^ +■■■ (3) 

Due to the symmetry only even states contribute (odd states are vanishing at 
a; = 0). On the other hand for the transition from a; = 1 to a:: = 1, we have 

(l|e-^^|l) = |(l|i?o)P e-^"^ + |(l|i?i)p e-^^^ + mE,)\^ e"^^^ + • • • (4) 

These transitions, Eqs. ©-lEJ, are sufficient to extract the low lying energy 
eigenvalues namely Eq , Ei and E2 ■ 

In the dilute gas approximation, a typical instantonic contribution to the 
transition in Eq. ^ is found to be composed of i + 1 instantons and i anti- 
instanton as shown in Fig. |21 which we call, for later use, M^'^. It should be 



x(t) 

+1 



\ 



ti 



\r^ 




Figure 2: The multi instanton contribution to the transition from x = to x = 1. The 
solutions with positive slope are called instantons whereas those of negative are called anti- 
instanton 

stressed that this contribution has 2* different configurations which should be 
taken into consideration during calculations. 

For the case of the transition, Eq. Q), one finds that it is composed of 
i— instantons and i anti-instantons as shown in Fig. O which we call Mf^^"^. 
This contribution has 2*~^ different configurations. Also there is a contribution 
coming from the trivial solution {x{t) = 1) that should be taken into consider- 
ation. 

It is worthy to note the following remarks. Firstly, in evaluating the transi- 
tion drawn in Figs. H2I3|I we have to compute the fluctuations over these paths 
consisting of strings of instantons-anti- instantons. However, because instantons 
are quite localized (with extension of order uj~^), we can neglect the fluctuations 
around their bodies. In other words, the quantum fluctuations would be effec- 
tive only around the straight line paths in Figs. H2I3|I . Along each step of the 
straight line paths, the particle resides at the minima of the potential (maxima 
for the inverted one) shown in Fig. ^ These fluctuations can be treated like the 
case of the harmonic oscillator. The crucial point here is to take care that the 
frequencies corresponding to each step may be different and can not be taken to 
be equal. Secondly, to compute the transition probability amplitudes given by 
Eqs. ©-0 one should sum over all the possible multi-instanton-anti- instanton 
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Figure 3: The multi instanton contribution to the transition from x = 1 to x = 1. The 
solutions with positive slope are called instantons whereas those of negative are called anti- 
instanton 



configurations. An elegant method for avoiding this comphcated sum has been 
developed and implemented in Rcfs. [HIIZ] to the case of quantum mechanical 
system with degenerate classical minima (exemplified by the symmetric double 
well) giving correct results for the low lying energy spectra and also non trivial 
information about the wave functions. The method in Refs. jS] [7] is based on 
saturating the transition probability amplitudes from the minimum to itself by 
the trivial solution, while for the transitions between two successive minima by 
one-instanton. However to extract energy eigenvalues one should have a priori 
some knowledge about the pattern of the correct energy spectra. Since in our 
present work we do not have such a priori knowledge about the correct pattern 
for the energy spectra for the symmetric triple well with non-cquivalcnt vacua, 
thus the sum over the multi-instanton-anti-instanton configurations can not be 
avoided. 

Using a technique similar to that used for the double well with non-cquivalcnt 
vacua [5], we rearrange the n- instanton integrals in such a form that recursive 
relations can be given for all these integrals. With the help of these relations the 
instanton sum can be easily managed. As stated before, only two transitions 
are required to be calculated, namely that between the minima at a; = and 
X = 1 and the other between a; = 1 to a; = 1. 

The contribution of the Mf has the integral: 
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where {N — -v/^) is the normalization constant, A is the exponential of the 
classical Euclidean action for one instanton and S = "'-""" . The factor K 
is calculated by matching the one instanton contribution. The one instanton 
contribution for the potential ,Eq. Q, can be found in Refs. 00], but with 



paying attention to used different conventions. 

Now tfiis integrai can be put in an cquivaicnt form: 
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wfiere B = V2K A . 

For the transition from x = 1 to a; = 1, the corresponding integral is 
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which can also be put in an another equivalent form 
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Carrying out the sum over the multi-instantons is rather involved, but can 
be simplified by studying the basic integral of the form 
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,,,.(£+^(£z^. 



(9) 



Integrating by parts we get the following recursive relations (here we drop the 
common factor e~2 2 "^ in the definition of I{n, m)) 
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With the help of these recursive relations, Eqs. pOfl . and taking into account the 
correct counting together with the starting initial integral /(O, 0) = -j e^ ~ — e^^ ~ 
we get (here we keep the common factor in the definition of /(n, m)) 

g-i^oTV- / TO + n-i 
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where (™) are the binomial coefRcients. 

For the transition from the minimum at a; = to the one at a; = 1 we need to 
evaluate the sum J^'ilo ^(*' *) '^hich is still hard to perform using Eq. pi|) . To 
facilitate evaluating the sum, we define S(n, m) — X]i=o -^(^^ + i,m + i) and let 
Sj{n,m) be the coefficient of S{n,m) associated with the terms - — y~ e^T. 
Using the recursion relations, Eqs. H1U|I . one gets 
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= S'iVi(n+ l,m), 

= S^{i, m + i — n) where n < i, 
= S^ (n + i — m,i) where m < i. 



(12) 



By letting a^ = S^{i, i), then with the aid of the recursion relations, Eas. p2|) . 
we get 

«m=«tiT;^a±, (13) 

consequently the sum over the multi-instantons becomes 
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The coefficient af corresponds to the coefficient of linear combination of two 
exponential (note that the factorial factor is removed) namely that 
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where a± are found to be 
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The coefficients C+ and C- can be then determined from 
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Thus from Eq. 111|) we get 
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leading to C+ = — -, — i ^ and C_ = 0. A similar treatment can be done for 

2\/l + (ro) 

a~ leading to the results C_ = ^ and C+ = 0. Using these results 

together with Eq. p4|) . one finally obtains for the odd multi-instantons sum: 

E Mfdd = /' [e-^o T _ ,-s. T] ^ (20) 
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with 
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Now the transition probability amplitude between the minimum at x = 1 to 
itself, after including the contribution of the trivial solution (x(i) = 1), becomes 



(l|e-^^|l) = TVe'^ +E Mf^en = ^,-^ + ^ ^ j(, _ 1,,), (22) 

which can be evaluated by the same procedure applied for the transition between 
the minimum at x = to that at a; = 1. For this purpose we let A^ = Sif{i, i+1) 
where Sf{i,i + 1) defined as before. Using the recursion relations in Eqs. H12|l 
we obtain 

A+ = f [Ati-A++i], (23) 

The coefficient A^ corresponds to the coefficient of linear combination of two 
exponential (note that the factorial factor is removed) namely that 

D+e^+ +D^e^-, (24) 



where X± are found to be 
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Again using Eq. (fTT|l we get 

and 
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which from Eqs. 12()ll lead to D+ = ^ — -i — , ^ and D_ = 0. A similar 

treatment can be done for A^ = S^_i{i — 1, i) leading to the recursive relation 
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After a lengthy algebra we finally get the instantonic contribution 
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where Eq and E2 as before, while Ei — ^. 

To sum up wc finally get our formulae for the energy spectra for the triple 
well, arranged in ascending order from the ground state,: 
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where B for the potential in Eq. (^ is: ^== w'^/^ e~^. 

Our formulae, Eqs. (|31|l . are completely different from that of Refs. |310]. 
In the limit of large w, the spectra tend to the limit Eq ^' ^, E2 ^ Ei = cu, 
i.e. a singlet and a degenerate doublet. 

To confirm our formulae we resort to numerical solution of Schrodinger equa- 
tion for the potential Eq. |^. In this letter we present the results for the energy 
differences in Tableland for the individual energy eigenvalues in Table El The 
details of the numerical method can be found in Ref. "^ . Moreover our results 
are consistent with what one could predict based on simple quantum mechanical 
approach to the problem Ref. JOI ■ 
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3.552713679 x lO^^^ 

2.135638335 x lO^^^ 



Table 1: The numerically calculated energy differences (AS""™) against those (AE^^^) 
predicted using instantonic approach. Notice that AEij = Ei — Ej . 
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Eq 14.178009913879056439 
El 27.856799752355288816 
E2 27.861522781409804929 



24.211602974742920912 
47.991273870020627020 
47.991274780080902476 



34.22366585434764411968736 
68.03469230180201792058325 
68.03469230190387712184547 



uj = 9Q 



L0= 110 



Eq 44.229945163238886600375384901 
El 88.056728233456712807586843449 
E2 88.056728233456721757979271787 



54.233802436877662843887532359 
108.070096366408926121501688376 
108.070096366408926122186175835 



Table 2: The numerically calculated energy eigenvalues. Numerical results, for largo u), show 
singlet and doublet structures [Eq — > ^, E2 — > -Ei = u)). 



In this letter we have emphasized the power of the application of the instan- 
ton method, when a proper treatment is used, and shown that the dilute gas 
approximation is sufficient. We have confirmed our splitting energy formulae 
numerically. As a final comment we expect that for a potential well with non- 
equivalent vacua, the multiplicity of the splitting is related to the number of 
equivalent vacua connected by the discrete symmetry of the potential. 
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